AD  682209 


ARPA  ORDER  NO.  180-1 


m 

MEMORANDUM 
RM- 5710- ARPA 

DECEMBER  1968 


FAR-FIELD  INTENSITY  DISTRIBUTION 
FROM  A  DIFFRACTING  APERTURE 
IN  A  TURBULENT  MEDIUM 


H.  T.  Yura 

•"D  D  r* 


rv 


7 


|1(  FEB  131968 

Ui/iSESED  U 


PREPARED  FOR 

ADVANCE  RESEARCH  PROJECTS  AGENCY 


1U 


R-pnD 


SANTA  MONICA  '  CALIFORNIA* 


fteprodocfd  by  tht 

CLEARINGHOUSE 

(■-  r  fedoral  Scientific  &  Technical 
Information  Springfield  Va  22151 


MISSING  PAGE 
NUMBERS  ARE  BLANK 
AND  WERE  NOT 
FILMED 


rAM c  m 


I 


MEMORANDUM 

RM-5710-ARPA 

DECEMBER  1968 


PAR-FIELD  INTENSITY  DISTRIBUTION 
FROM  A  DIFFRACTING  APERTURE 
IN  A  TURBULENT  MEDIUM 

H.  T.  Yura 


Tin-  rc'eari'h  i*  supported  In’  tlic  Advanced  llc»oareh  Projects  Aficney  under  Contract 
Nii.  I)  AIK  15  67  C  01  II.  View*  or  conclusion'  eoiil.-iined  in  this  study  should  r.ol  lie 
interpreted  ;is  representin':  the  official  opinion  or  policy  of  AIII’A. 


1)1  STH I  Bl’TION  STATEMENT 

This  document  has  horn  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 


p-fl  n  D 


I  »  O  *  N  I  * 


*0««* 


This  study  is  prosmlrd  a*  a  ruuipctrnl  Irralmnit  of  I h«*  snlijn  t.  worthy  of  piili- 
liratiou.  The  Hand  Corporation  row-ln1*  for  tin-  tpiality  of  the  rrsrnrrli.  without 
iKTossarilv  «Miilor»in"  the  opinion-  and  rou<-lii“ion*  of  tin*  authors. 


Published  by  The  RAND  Corporation 


lii 


PREFACE 


This  Memorandum,  prepared  for  the  Advanced  Research  Projects 
Agency,  is  part  of  a  study  of  those  phenomena  which  affect  the  per¬ 
formance  of  optical  or  infrared  reconnaissance  and  guidance  equipment. 
The  objective  of  these  studies  is  to  provide  sufficient  understanding 
for  the  system  analyst  to  compute  performance  estimates  under  various 
operational  conditions. 

A  quantitative  understanding  of  the  effect  of  atmospheric  turbu¬ 
lence  on  a  beam  of  light  of  finite  cross  section  is  required  for  the 
prediction  of  the  performance  of  various  devices  employing  lasers  for 
target  acquisition  or  guidance  in  tactical  missions.  Such  applications 
are  characterized  by  near-horizontal  propagation  paths  near  the  ground 
of  the  order  of  one  to  tens  of  kilometers  in  length.  This  Memorandum 
obtains  the  modified  intensity  distribution  in  the  far  field  of  a  fi¬ 
nite  transmitting  aperture  due  to  a  turbulent  medium.  These  results 
should  be  of  use  to  those  interested  in  tactical  applications  of  laser 
range  finders,  laser  line  scanners,  and  the  various  guidance  systems 
employing  an  illuminating  beam. 
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SUMMARY 

The  far-field  Intensity  distribution  at  optical  frequencies  of  an 
initially  plane  wave  from  a  finite,  circular,  source  aperture  is  ob¬ 
tained  as  a  function  of  range  and  angle  for  various  values  of  the  index 
structure  constant.  Describing  the  turbulence-induced  index  of  refraction 
fluctuations  by  the  Kolmogorov  spectrum,  it  is  found  that  the  time 
average  of  the  peak  radiant  intensity  in  the  Fraunhofer  region  of  a 
finite  source  aperture  decreases  with  range  at  a  much  faster  rate  than 
the  intensity  calculated  from  absorption  and  scattering  by  molecules  and 
particulate  matter.  An  expression  is  derived  for  this  intensity  as  the 
sum  of  two  terms:  one  which  represents  the  exponential  decay  of  the 
energy  in  the  initially  coherent  beam;  and  the  other  which  represents 
the  complementary  growth  of  the  energy  in  an  incoherent  radiation  field, 
attributed  to  the  loss  of  coherence  by  scatterings  off  the  turbulence- 
induced  fluctuations  in  the  medium.  The  exponent  in  this  process  is 

shown  to  be  of  the  order  (R/R  (where  R  is  the  range  where  the 

c  C  -1 

average  part  of  the  field  is  down  by  a  factor  of  the  order  e  )  in 

contrast  to  the  usual  (R/Rv)  for  the  attenuation  by  molecules  and  dust 

(where  R^  is  the  visual  range).  The  beamwidth  at  half-power  is  shown 

to  increase  very  slowly  until  the  propagation  distance  reaches  a  range 

that  is  of  the  order  of  R  .  These  results  can  be  illustrated  by  the 

c 

application  to  a  beam  from  an  aperture  of  2-cm  diameter  at  a  wavelength 

of  0.6328}i.  For  moderate  daytime  turbulence,  R  ^  1.7  km  at  this 

c 

wavelength,  compared  to  ar.  R^  of  about  5  km  for  a  correspondingly 

moderate  visual  range.  For  small  values  of  l  (e. g. .  I  £  0.5  cm) 

o  o 

where  lQ  is  the  inner  scale  of  turbulence,  the  effects  of  turbulence 

dominate  over  most  of  the  range  of  interest.  For  this  case,  the  beam 

will  exhibit  some  structure  out  to  a  range  of  about  (1.0  -  1.65)  R 
.  .  c 

for  0.5  cm  c  1  0. 1  cm. 
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I.  INTRODUCTION 


In  the  presence  of  turbulence,  inhomogeneities  appear  in  the  me¬ 
dium  which  scatter  electromagnetic  waves.  The  electromagnetic  field 
diffracted  by  an  aperture  exhibits,  in  the  absence  of  turbulence,  a 
cl aracteristic  diffraction  pattern  with  well-defined  regions  of  maximum 
and  minimum  intensity.  The  purpose  of  this  Memorandum  is  to  calculate, 
at  optical  frequencies,  the  effects  of  the  turbulent  medium  on  the  dif¬ 
fraction  pattern.  It  is  shown  that  the  beam  pattern  flattens  and 
broadens  as  the  propagation  distance  increases.  Analytic  expressions 
are  derived  for  the  intensity  distribution  of  a  plane  wave  incident 
on  a  circular  aperture  as  a  function  of  propagation  distance  and  for 
the  parameters  that  characterize  the  turbulent  medium. 

The  results  obtained  here  have  application  in  predicting  the  per¬ 
formance  of  systems  employing  a  laser  for  the  detection  and  location  of 
targets.  The  assumption  is  made  that  the  parameters  that  characterize 
the  turbulent  medium  are  not  a  function  of  range.  This  is  a  good  as¬ 
sumption  for  horizontal  paths.  However,  these  results  may  also  be 
extended  directly  to  apply  to  conmunications  from  satellites,  etc. 

The  present  analysis  is  restricted  to  a  weakly  inhomogeneous  me¬ 
dium.  That  is,  it  Is  assumed  that  the  fluctuations  are  small  compared 
to  the  average  properties  of  the  medium.  Furthermore,  it  is  assumed 
that  the  characteristics  of  the  medium  do  not  change  appreciably  in  a 
period  of  oscillation  of  the  electromagnetic  field,  because  frequency 
spreading  (doppler  effects)  then  becomes  important.  At  near -infrared 
and  optical  wavelengths,  this  condition  is  satisfied  in  the  atmosphere. 
The  electromagnetic  field  under  consideration  has  a  time  dependence 

-ittit 

given  by  the  factor  e  .  In  this  case  the  time-dependent  wave  equa¬ 
tion  is  replaced  with  the  Helmholtz  equation  for  an  inhomogeneous  me¬ 
dium.  The  electrical  conductivity  and  magnetic  permeability  of  the 
medium  are  taken  to  be  zero  and  one  respectively.  Here,  only  the  case 
of  the  propagation  of  a  scalar  field  in  a  fluctuating  medium  is  dis¬ 
cussed.  Extension  to  vector  fields  is  straightforward.  It  is  noted 
that  the  results  obtained  here  apply  not  only  to  electromagnetic  fields, 


but  to  any  field  that  satisfies  the  wave  equation,  e.g. ,  acoustic 
waves  when  the  scale  of  turbulence  is  greater  than  the  wavelength. 
Furthermore,  it  is  assumed  that  the  diameter  D  of  the  aperture  is 
larger  than  the  inner  scale  of  turbulence  t .  Under  this  condition, 
effects  such  as  beam  steering  may  be  neglected. 

In  Section  II  the  physical  nature  of  the  electromagnetic  field  in 
a  turbulent  medium  is  reviewed.  A  previously  demonstrated  treatment 
of  the  field  as  a  sum  of  an  average  and  a  random  part  is  described. 

In  Section  III  the  time  average  of  the  electromagnetic  field  is 
discussed  and  applied  to  the  field  diffracted  by  an  aperture.  In 
Section  IV  the  random  field  is  discussed,  and  in  Section  V  the  in¬ 
tensity  distribution  as  a  function  of  range  and  field  angle  is 
obtained.  It  should  be  noted  that  these  results  are  derived  on  the 
assumption  of  the  single  scattering  approximation  and  therefore  may 
not  be  accurate  for  R  »  R£  (R£  is  the  range  where  the  average  part 
of  the  field  is  down  by  a  factor  of  the  order  e  *).  Curves  depicting 
the  beam  pattern  (based  on  the  Kolmogorov  spectrum  for  the  index  of 
refraction  fluctuations)  as  a  function  of  various  parameters  are  given 
It  is  shown  that  the  beam  pattern  flattens  and  broadens  as  the  propa¬ 
gation  distance  R  Increases.  In  particular,  when  R  >  Rj^  (defined 
below)  the  beam  pattern  is  essentially  flat. 

The  analytic  results  obtained  are  valid  for  propagation  distances 
2 

R  >  Rj  ■  D  /X,  where  D  and  X  are  the  diameter  of  the  aperture  and  wave 
length  of  the  field  respectively;  i.e.,  the  results  are  valid  in  the 
Fraunhofer  region.  It  is  assjmed  here  that  R.  <  R  •  In  this  case, 
the  modification  of  the  nerr-field  diffraction  pattern  (Fresnel 
region)  is  entirely  negligible. 
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II.  GENERAL  CONSIDERATIONS 


In  this  section  the  qualitative  nature  of  an  electromagnetic  wave 

propagating  In  a  weakly  Inhomogeneous  medium  Is  reviewed  briefly. ^ 

U,  a  typical  component  of  the  field,  can  be  represented  as  the  sum  of 

• — 

two  terms:  U  (R) ,  the  average  field  defined  by  U  ■  U:  and  a  ran- 

dom  field  U  where  Uran  ■  0.  As  the  field  progresses  through  the 

medium,  energy  Is  transferred  from  the  average  part  to  the  random  part 

of  the  field  so  that  energy  Is  conserved.  For  propagation  distances  R 

much  less  than  a  critical  distance  R  (defined  below)  the  field  Is 

c 

primarily  unscattered;  the  scattered  component  Is  of  the  order  n^  where 

n^  Is  the  fluctuating  part  of  the  Index  of  refraction  (l.e.,  n  ■  1  +  n^, 

|n^|  <<  1,  n^  ■  0).  For  propagation  distances  much  greater  than  R£ 

the  field  Is  essentially  random  since  the  average  field  Is  exponentially 

small  (for  R  >>  R^  the  field  has  undergone  multiple  scatterings  off  the 

random  Inhomogeneltles  of  the  medium  and  Is  completely  random) . 

For  R  <<  R  ,  the  total  field 
c 


U(ft> 


U  (ft) 

avg' 


U  (ft) 
ran' 


(1) 


where,  to  terms  up  to  second  order  in  n^,  Uflv  (R)  is  formally  the  same 
as  the  field  in  the  absence  of  turbulence  (l.e.,  for  n^  =  0),  except 
that  the  propagation  wave  vector  k  Is  not  given  by  Its  value  In  the 
absence  of  turbulence  (kQ)  but  is  given  by^ 


k 


k 


o 


Bn(0) 

2 


k2 

+  2? 


>00 

n  o 


(2) 


where  B  (0)  Is  the  correlation  function  of  the  Index  of  refraction  fluc- 
n  ** 
tuatlon  evaluated  at  r  ■  0,  l.e.. 


A  bar  over  a  quantity  signifies  the  ensemble  average  of  the  quan¬ 
tity. 

** 

It  Is  assumed  here  that  the  medium  is  statistically  homogeneous 
and  Isotropic. 


-  Jf  • 


Bn(r)  -  n1(r1)n1(r2) 


Bn(rl  “  r2} 


Bn(r)  where  r  -  ri  “  r2^ 


/o 

dr  —  B(r)e 
r  n 


ik  •  r 
o 


It  can  be  seen  that  $  (k  )  (and  hence  k)  has  a  positive  imaginary 
)  no 

part. 

The  random  part  of  the  field,  to  lowest  order  in  n^,  is  given  by 
the  first  Born  approximation: 


U  (R)  -  U  (R)ik(R) 
ran  o  rl 


where  U  (R)  is  the  field  in  the  absence  of  turbulence,  and  ip,  (R)  is 

°(2)  1 
given  by 


fc- 

,  (R)  -  - /n  (R')U  (R') 

*  n_»T  /n\  J  *  ° 


2ttU  (R) 
o 


ik  Ir-3' 

o' 

I 

|$  -  R' 


where  the  integration  volume  is  that  region  of  space  where  n^R)  i  0 
(the  aperture  being  located  at  R  ■  0). 
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III.  THE  AVERAGE  FIELD 


Consider  an  expanding  spherical  wave.  In  the  absence  of  turbulence 
this  field  Is  given  by 


ik  R 

Uo(R)  "  c  ^TT-  (6) 

where  C  Is  an  arbitrary  constant.  For  R  <<  R  the  average  part  of  the 
spherical  wave  Is  given  by 


lkR 


(7) 


Let 


*  ■  i(k  -  k  )R 
o 


(8) 


Then 


U 

avg 


<R) 


ik  R  r  ^ 

Ce _ Re<|>  +  ilm$  I 

R  l*  J 


(9) 


To  lowest  order  in 


Re  <J> 


k3R 

o 

2x 


Imf  (k  ) 
n  o 


(10) 


and 


6 


Im$ 


rB(o>  ] 

[—  +  T  lm*n(ko>JkoR 


(ID 


Thus,  In  Che  presence  of  the  random  inhomogeneities,  the  average  compo¬ 
nent  of  the  field  possesses  a  positive  attenuation  coefficient.  The 
quantitative  results  given  by  Eqs.  (2)  and  (5)  are  obtained  by  a  pertur¬ 


bation  method  valid  for  R  <  R 


The  critical  distance  R  is  obtained 

c 


from  the  requirement  that  the  relative  magnitude  of  the  correction  to 
the  field  in  the  absence  of  turbulence  be  small.  It  can  be  seen  that 
this  requirement  leads  to  the  condition  that^ 


| Re$ |  <<  1 


(12) 


The  critical  distance  R  is  determined  by 

c  1 


|Re$(Rc)|  -  1 


(13) 


Thus,  R  is  the  propagation  distance  where  the  average  component  of  the 

C  -1 
field  is  down  by  a  factor  of  the  order  e 

It  follows  Immediately  from  the  optical  theorem^  (or  from  energy 

conservation)  that 


-  2Re$(R)  -  |  d'1(R)  | 


Re^(R)J2  +  |lm<fl(R)J: 


(14) 


where  is  given  by  Eq.  (5).  Hence  R£  may  also  be  determined  from 


1/2  i*1(Rc)r 


i 


(15) 
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1/2 

Furthermore,  if  (X  R)  >>  £  ,  where  i  is  the  inner  scale  of  turbu- 
*  o  o  o 

lence,  it  can  be  shown  that(2) 


(Reit<1)2~  (Imi^)2 

In  this  case  R£  may  be  determined  from 


(16) 


[Re^(Rc)]2  -  1,  (\>R)1/2  »  *0 


(17) 


Choosing  the  Kolmogorov  spectrum  to  represent  the  index  of  refrac¬ 
tion  fluctuations,  it  can  be  shown  that  for  a  spherical  wave  and 

(X  R)1/2  »  i 
o  o 


( Re i|).  ) 2  a:  0.13  C2k7/6  RU/6 
1  no 


(18) 


** 


2  (2) 
where  C  is  the  index  structure  function.  Strictly  speaking,  Eq. 

n  2 
(18)  is  valid  for  propagation  paths  where  C  is  independent  of  posi¬ 
tion  along  the  path,  e.g.,  horizontal  paths.  For  slant  paths  through 
the  atmosphere,  Eq.  (18)  (for  spherical  waves)  is  replaced  by 


(Rei^)2 


(18a) 


***  2 

where  h  ■  z  cos  a,  and  a  is  the  zenith  angle.  When  C  is  a  con- 

n 

stant,  Eq.  (18a)  is  identical  to  Eq.  (18).  In  this  paper  it  is  as- 
2 

sumed  that  C  is  a  constant.  Extension  to  the  more  general  case  is 
n 

straightforward. 


Typically,  0.1  cm  £  l  s  1  cm  for  horizontal  paths  near  sea  level. 
**  0 
See  Ref.  2,  Eq.  9.43. 

*** 

See  Ref.  2,  Chap.  8. 
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2 

The  values  of  C  typical  of  turbulence  in  the  first  few  hundred 
n  -1'  -2/3 

meters  of  the  atmosphere  are:  30  x  10  "cm  (strong  daytime  tur- 

-15  -2/3 

bulence  within  a  few  meters  of  the  ground)  3  '  10  cm  and 
-15  -2/3 

7  x  10  cm  (moderate  daytime  turbulence  «nd/or  strong  nighttime 

-15  -2/3 

turbulence),  and  0.3  x  10  cm  (very  weak  turbulence  which  occurs 

In  the  near  neutral  periods  at  dawn  or  dusk).  Figure  1  gives  R  as  a 

2  c 

function  of  wavelength  for  various  values  of  C  •  For  example,  for 
\Q  ■  0.6328u  (He-Ne  laser  wavelength)  It  is  found  that 


R  —  5.9  km, 
c  * 

^  3.1  km, 


i  1.7  km, 


—1.1  km, 


2!  0.48  km, 


r2  n  i  v  in-15  “2/3 

C  ■  0.3  x  10  cm 
n 


C2  -  1  x  10“15  cm“2/3 
n 


C2  -  3  x  10“15  cm  2/3 
n 


c2  -  7  x  I(f 15  cm_2/3 
n 


P2  in  „  in_15  -2/3 

C  »  30  x  10  cm 
n 


(19) 


For  plane  waves,  Rc  is  equal  to  0.61  times  the  values  for  a  spherical 
wave  given  by  Eq.  (19). 

Next,  consider  the  average  field  diffracted  by  an  aperture.  In 

the  Fraunhofer  region,  the  field  in  the  absence  of  turbulence  is  given 
(3) 

by  (it  is  assumed  here  that  Aq/D  <<  1  where  D  is  a  characteristic 
length  of  the  aperture) 


(20) 


where  u  is  the  magnitude  of  the  field  (assumed  constant)  in  the  aper- 
turn,  p  is  the  vector  from  the  origin  to  an  arbitrary  point  in  the 
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Fig.  I— The  critical  distance  Rc  as  a  function  of  wavelength 
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aperture,  k'  Is  the  wave  vector  of  the  diffracted  wave  (differing  in  di 
rection  only  from  kQ  which  is  assumed  nearly  normal  to  the  aperture) , 
and  dA^  is  the  projection  of  an  area  element  of  the  aperture  normal  to 
the  direction  of  the  initial  (undiffracted)  light  wave. 

In  the  presence  of  the  random  medium  the  average  field  is  given 
by  Eq.  (20),  with  kQ  replaced  by  k  where  k  is  given  by  Eq.  (2).  For 
a  circular  aperture  of  radius  a,  from  Eq.  (20) 


U 

avg 


(R.O) 


(21) 


where 


e  •  ka  \ 

X  -  R6  |  (22) 

sin  0  -  0  f 


and  k  is  given  by  Eq.  (2). 
To  lowest  order  in  n^ 


Jv  I 

(“  <>') 

/Ji(ax0>\ 

)  r  i 

j„xp  -•i|*1<«>|2 

i 1  1 

\  *  ) 

V  x„  ) 

+  O  n?  I  +. 


(23) 


where  xq  ■  kQ0  and  i/^(R)  is  given  by  Eq.  (5).  The  terms  in  the  braces 
in  Eq.  (23)  represent  the  diffracted  field  in  the  absence  of  turbulence, 
and  the  effects  of  the  random  medium  are  represented  by  the  exponential 
term.  The  analytic  expression  for  given  by  Eq.  (5)  is  valid  for 
R  <  Rc«  On  physical  grounds  it  is  expected  that  the  qualitative  nature 
of  the  average  field  remains  the  same  for  R>  R£;  i.e.,  it  is  damped. 

In  particular,  for  large  propagation  distances,  Ufl  »  0,  all  of  the 
field  energy  having  been  transfered  to  the  random  field. 
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IV.  THE  RANDOM  FIELD 

The  random  field  for  R  <<  R^  is  given  by  Eq .  (5).  In  this  expres¬ 
sion  the  quantity  (Jq(R)  Is  the  field  that  is  present  in  the  absence  of 
turbulence.  In  the  case  of  Fraunhofer  diffraction  by  an  aperture  (lo¬ 
cated  at  the  origin  of  coordinates),  the  random  field  <Jj j  is  thus  given 
by  Eq.  (5)  with  Uq  given  by  Eq .  (20).  In  the  integral  of  Eq.  (5)  it 
is  seen  that  an  integration  of  Uo(R')  over  all  values  of  R'  is  required. 

Equation  (20)  gives  the  diffracted  field  in  the  Fraunhofer  region 
2  2 

(R  >  D  /X).  For  R  <<  D  /X,  the  Fresnel  expression  for  the  diffracted 

field  should  be  used.  In  this  report  interest  is  in  large  propagation 
2 

distances  (R  »  D  /X);  hence  negligible  error  is  introduced  in  the  vol¬ 
ume  integration  of  Eq .  (5)  by  using  Eq.  (20)  for  all  values  of  R. 

The  integral  obtained  by  substituting  Eq .  (20)  into  Eq.  (5)  pre¬ 
sents  great  mathematical  difficulty.  However,  if  we  consider  the  case 

k^KR  «  1  (24) 

(where  K  is  the  wave  number  of  the  power  spectrum  of  the  index  of  re¬ 
fraction  fluctuation)  it  can  be  shown  that  the  resulting  integral  ex¬ 
pression  obtained  reduces  to  that  given  by  Tatarski  for  a  spherical 

*  1/2 

wave.  Furthermore,  long  propagation  distances  where  (X  R)  »  l 

o  o 

are  of  interest.  In  this  case  the  spectrum  of  the  correlation  function 
for  the  incoherent  field  fluctuation  is  concentrated  near 

K  2*/(XqR)1/2  (25) 

(i.e.,  the  refractive  index  inhomogeneities  with  scales  of  the  order 
1/2 

(X  R)  make  the  largest  contribution  to  the  field  fluctuations  of 
°  ** 

the  wave).  From  Eqs.  (24)  and  (25)  we  find  that  in  the  case  of 

* 

See  Ref.  2,  Chap.  9. 

** 

See  Ref.  2,  Chap.  7. 
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R  »  Rx  -  2-  (26) 

Ao 

Tatarskl's  results  for  a  spherical  wave  may  be  used  for  the  integral 

-*  2 

expression  of  f^(R)  •  For  a  wide  range  of  parameters,  R,  «  R£.  For 

X  ■  0.6328y,  and  a  ■  1  cm,  R.  —  600  m,  while  R  >  1  km  (see  Eq.  (19)). 
o  i  c  2 

Results  derived  here  apply  to  the  case  R  >  R^  ■  D  /\q.  For  R  < 

a  more  thorough  analysis  of  the  integral  expression  for  ^  (i.e.,  Eq. 

(5))  with  Uq  given  by  Eq.  (20)  is  needed.  No  qualitative  calculations 

for  R  <  Rj  are  attempted  here. 
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V.  BEAM  PATTERN 

In  this  section  the  results  of  Sections  III  and  IV  are  used  to 
obtain  the  diffraction  pattern  of  a  circular  aperture  in  the  presence 
of  the  inhomogeneous  medium.  First,  an  approximate  expression  for  the 
diffracted  power  valid  for  R  ~  Rc  is  derived.  To  do  this  it  is  as¬ 
sumed  that  the  average  radiated  power  can  be  written  as  a  sum  of  two 

terms  for  all  values  of  R:  an  average  part  (proportional  to  ju  |^), 

I  1 2 

and  a  random  part  (proportional  to  |U  |  ).  This,  in  effect,  assumes 
that  the  phase  of  the  random  part  of  the  field  is  not  correlated  to 
the  phase  of  U 


_avg 


This  is  true  both  for  R  »  R  and  for  R  «  R  : 

c  c 


i.e. ,  U  U  "0  for  both  of  these  limiting  cases,  and  it  is  also 
ran 

assumed  to  be  true  for  R  ~  Rc>  In  order  to  obtain  the  angular 

dependence  of  the  beam  pattern  it  is  noted  that  the  angle  of  scattering 

of  the  field  by  the  refractive  index  inhomogenieties  is  of  the  order 

\/jf  ,  where  £  is  the  inner  scale  of  turbulence  (see  Ref.  2).  The 

assumption  that  the  scattering  angle  is  of  the  order  X/£q  i8  consistent 

with  the  observation  that  stellar  images  even  under  poor  seeing  con- 

(4) 

ditions  exhibit  a  finite  blur  circle.  In  the  absence  of  turbulence 
almost  all  of  the  diffracted  power  is  contained  in  a  cone  of  the  order 
X/D.  Since  £q  <  D  (typically,  0.1  cm  5  M  1  cm),  it  is  seen  that  most 
of  the  scattered  power  will  be  contained  in  a  cone  with  angular  aperture 
of  the  order  X/£q.  Hence,  to  obtain  the  beam  pattern  we  assume  that 
the  scattered  radiation  is  contained  in  this  cone.  This  neglects  the 
relatively  small  amount  of  radiation  contained  in  the  sidelobes.  In 
view  of  other  approximations  involved  in  the  calculation  and  uncertainties 
in  the  underlying  statistical  description  of  the  index  of  refraction 
fluctuations,  it  is  felt  that  this  approximation  is  justifiable. 

Therefore,  let  (dT<  «  2it  sin  6  d0) 


dP 


dP  +  dP 
avg  ran 


(27) 


where 


dP 


V 


avg 


—■  A^(Xa)  e‘£(R)dr> 


(28) 
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P_  Is  the  power  In  the  aperture,  and  A-(x)[*  2J-(x)/x]  is  the  Lotnr.el 

U  (5)  1  1 

function  of  the  first  kind  whose  square  gives  the  beam  pattern  in 
the  absence  of  turbulence,  and 


dP 


ran 


C_ 

2 n 


(1  -  e 


"£<R))dn 


(29) 


The  constant  C  in  Eq.  (29)  is  determined  by  requiring  that^  J  dP  ■  Pq 
(the  power  in  the  aperture),  which  yields  that  C  ■  2Pq  (?)  •  - 

quantity  0R  is  determined  by  the  requirement  that  for  small  R,  the 
power  given  by  Eq.  (27)  gives  the  perturbation  result.  Doing  this  we 
obtain 


Hence, 


f(R)  -  2[Re*lkR)]2 

■  0.26cV/6R11/6  (for  the  Kolmogorov  spectrum) 
from  Eqs.  (27) -(30) 


(30) 


where 


(31) 


For  R  »  Rr ,  dP  ■  Pgdfl}  the  average  field  being  exponentially 

small,  the  quantity  in  the  braces  in  Eq.  (31)  represents  the  modified 

beam  pattern  due  to  the  effects  of  turbulence.  Denoting  this  quantity 
2 

by  A1(eO)  F(R,0),  we  have 


A‘(e6)  F(R,6) 


A^(e0) 


c  *  n  a 


2  2 
V  A{(c0) 


(32) 
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2 

where,  to  lowest  order  in  n^,  e  ■  kQa  ■  2na/X  «  1.  The  quantity  in 
the  brackets  in  Eq.  (32)  thus  represents  the  effects  of  the  turbulent 
medium  on  the  beam  pattern. 

A  quantity  of  interest  is  the  dependence  of  the  radiant  intensity 

2 

on  range  for  9  ■  0.  Figures  2  and  3  are  a  plot  of  A^(0)  F(R,0)  as  a 
function  of  propagation  distance  R  for  the  Kolmogorov  spectrum  for  var¬ 
ious  values  of  the  index  structure  constant,  £  and  X  -  0.6328p,  and 
1.06|i.  In  addition  to  turbulent  scattering,  there  is  an  exponential 
loss  of  intensity  due  to  both  absorption  and  scattering  by  molecular 
and  particulate  constituents  of  the  atmosphere.  This  attenuation  is 
given  by  exp[-R/Rv]  where  R^  is  the  visual  range. 

Commonly  occurring  values  of  Ry  are  1-5  km.  Included  in  Figs. 

2  and  3  is  a  plot  of  this  attenuation  for  R^  *■  1  and  5  km.  It  is  seen 
that  even  with  moderate  turbulence  and  D  *  2  cm,  the  attenuation  of 
the  beam  will  be  dominated  by  turbulent  scattering  when  £q  £o.5  cm. 

For  larger  values  of  D  within  the  range  of  interest  it  is  expected 
that  turbulent  scattering  will  be  dominant  for  £Q  ~  1  cm*  In  Figs. 

4  and  5  the  pattern  (i.e.,  Eq.  (32))  is  plotted  as  a  function  of 

t(-  €0)  for  various  values  of  propagation  distance  R,  t  ,  and 

2  -15  -2/3  ° 

Cno  ■  3  x  10  cm  for  a  2-cm  diameter  aperture.  The  quantity 

2 

A  (t)  gives  the  beam  pattern  in  the  absence  of  turbulence.  The  beam 
1  2 

pattern  for  other  values  of  C'  (for  fixed  l  ,  X,  and  D)  may  be 

n  o  226/11 

obtained  from  Figs.  4  and  5  by  replacing  R  by  R[C  / C  j  '  where 

2  -15  -2/3  no  n 

C  ■  3  X  10  cm  .  The  sharp  discontinuities  shown  in  Figs.  4c 
no 

and  5c  result  from  the  assumption  that  the  scattered  power  is  contained 
in  a  cone  of  half-angle  equal  to  X/£q  and  are  not  physically  significant. 

Another  quantity  of  interest  is  the  half-power  width.  That  is, 
that  value  of  9  where  the  radiant  intensity  is  one-half  its  value  at 
0  ■  0.  From  Eq.  (32)  the  angle  0^^  at  half-maximum  is  given  by  the 
solution  to 


Note  that  A^(0)  ■  1;  in  the  absence  of  turbulence  the  normalized 
Fraunhofer  beam  pattern  is  independent  of  range. 


Relative  radiant  intensity  (9  s  0) 
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I  0 


0.1 


0.01 


0.1  10  10  100 

Range  (  Km ) 


Fig.  2— The  dependence  of  radiant  intensity  on  range  for 
0=0  (D  =  2  cm,  X=  0.6328 fi)  and  various  values  of  JL0 


Relative  radiant  intensity  (9  =  0) 
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Range  (Km  ) 


Ronge  (Km) 


Range  (Km) 


Fig.  3— The  dependence  of  radiant  intensity  on  range  for 
9=  0  ( D  =  2  cm,  X=  1. 06 fj.)  and  various  values  of  /0 
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Rearranging  terms  in  Eq.  (33) , 


Equation  (35)  has  been  solved  numerically  for  t^  38  3  function  of 
range  and  is  plotted  in  Figs.  6  and  7  for  various  values  of  C^, 
and  X.  As  R  inr  eases,  t^2  increases  slowly  until  the  propagation 
distance  reaches  a  range  that  is  of  the  order  of  R£ ,  where  it  increases 
rapidly  to  a  limiting  value  equal  to  TiD/ t^.  The  curves  in  Figs.  6  and 
7  are  for  range  values  such  that  t^2  <  "D /lQ,  the  quantity  tQ  1.62) 
is  the  half-power  width  in  the  absence  of  turbulence. 

The  diffracted  beam  flattens  and  broadens  as  R  Increases.  That 
is,  for  large  enough  R,  the  beam  pattern  is  contained  in  a  cone  angle 
of  the  order  X/ For  values  of  t(=  60)  near  the  central  maximum  of 
the  beam  pattern  (~  tt)  ,  R^  is  defined  as  that  range  such  that  for 
R  >  R^  the  radiant  intensity  never  varies  by  more  than  one-half  its 
maximum  value.  For  R  >  R^  the  effect  of  the  transmitting  aperture 
is  lost.  We  note  that  t^^2  »  1  at  this  range;  hence,  from  Eq.  (35) 
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Fig.  6— The  normalized  half-power  width  t1/2/t0  as  a  function 
of  range  for  X  =  0.6328/2.,  D  =  2  cm,  various  values  of 
the  index  structure  constant,  and  i0 


to 


Fig.  7 — The  normalized  half-power  width  t,/2/t0  as  a  function 
of  range  for  X=  1.06/x,  D  =  2  cm,  various  values  of  the 
index  structure  constant,  and  J0 
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Thus,  as  R  increases,  t^^  increases  from  its  value  in  the  absence  of 
turbulence  (tQ  =  1.62),  tending  rapidly  to  its  limiting  value  (ttD/jIo) 
as  R  -*  R^.  Figure  8  is  a  plot  of 
as  a  function  of  wavelength  for  v. 
a  plot  of  the  coefficient  of  Rc  in  Eq.  (37).  This  factor  is  independent 
of  wavelength  and  turbulence  strength  and  is  only  a  function  of  l 
and  D. 

These  results  have  been  derived  assuming  the  single  scattering 
approximation.  It  is  expected  that  this  assumption  is  valid  for 
R  ^  Rc«  Therefore  the  asymptotic  behavior  (for  R  >  Rc)  shown  in  Figs. 

2  and  3  may  be  an  artifact  of  this  approximation. 

These  results  can  be  illustrated  by  the  application  to  a  beam 
from  an  aperture  of  2-cm  diameter  at  a  wavelength  of  0.6328a.  For 
moderate  daytime  turbulence,  Rc  =■  1.7  km  at  this  wavelength  compared 
to  an  R^  of  about  5  km  for  a  correspondingly  moderate  visual  range. 
However,  for  (  <^0.5  cm,  the  difference  in  functional  dependence  results 
in  a  crossover  range  with  the  effects  of  turbulence  dominant  over  most 


R^,  for  a  2-cm  diameter  aperture, 

rious  values  of  C^.  Figure  9  is 
n 


of  the  range  of  interest.  For  this  case,  the  beam  will  exhibit  some 

structure  out  to  a  range  of  about  (1.0  -  1.65)  R  for  0. 5  cm  ~  l  £  0. 1 

c  o 


cm. 


(Km)  Ru(km)  Ru  (km) 


2  A 


Fig.  8— The  maximum  range  RM  as  a  function  of  wavelength 
for  various  values  of  the  index  structure  constant, 

D  =  2  cm,  and  J0 


L 
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